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ABSTRACT 

Multicolor HST WFPC-2 images are used to study the central structure of the 
three Local Group galaxies M31, M32, and M33. PSF-deconvolution and modeling 
of image aliasing are required to recover accurate brightness profiles within r < 0"5. 
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The data present a study in contrasts that suggests different evolutionary histories. In 

M31, the nucleus is doublc-pcaked, as found by WFPC-1 and confirmed by WFPC-2. 
The dimmer peak, P2, is closely centered on the bulge isophotes to O'.'l, implying that 
it is the dynamical center of the galaxy. Directly on P2 lies a UV-bright compact 
source that was discovered by King ct al. (1995) at 1700 A. WFPC-2 images now 
show that this source is resolved, with ri/2 ~ 0.2 pc. It dominates the nucleus at 3000 
A, and its spectral energy distribution is consistent with late B-early A stars. This 
probable nuclear star cluster may consist of young stars and be an older version of the 
central cluster of hot stars that now sits at the center of the Milky Way, or it may 
consist of heavier stars built up from collisions in a possible cold disk of stars orbiting 
P2. Aspects of its formation remain highly problematic. In M32, new images show 
that the central cusp continues to rise into the HST limit with slope 7 0.5, and the 
central density po > lO^M© pc~^. The V — I and U — V color profiles are essentially 
flat, and there is no sign of an inner disk, dust, or any other structure. This total lack 
of features seems at variance with a nominal stellar collision time of 2 x 10^*^ yr, which 
implies that a significant fraction of the light in the central pixel should come from 
blue stragglers. This discrepancy is eased but not completely removed if the stellar 
population is young (2 x 10^ yr). The stubborn normalcy of M32 at tiny radii may 
be emerging as an important puzzle. In M33, the nucleus has an extremely steep 
7 = 1.49 power-law profile for 0''05 < r < 0"2 that appears to become somewhat 
shallower as the HST resolution limit is approached. The profile for r < 0^'04 can 
be described as having either a 7 « 0.8 cusp or a small core with Vc ~ 0.13 pc. The 
central density is po > 2 x 10^ M0 pc~^, and the implied central relaxation time is 
only ~ 3 X 10^ yr, indicating that the nucleus is highly relaxed. The accompanying 
short collision time of 7 x 10^ yr predicts a central blue straggler component that is 
quantitatively consistent with the strong V — I and B — R color gradients seen with 
HST and from the ground. When combined with the Galaxy, the nuclei of the Local 
Group show surprisingly similar radial density profiles but divide into two classes 
according to velocity dispersion and black hole content: M31, M32, and the Milky 
Way are dominated dynamically (and stabilized against relaxation) by massive central 
black holes, while M33 lacks a dominant black hole. An obvious hypothesis is that a 
sizable stellar spheroid (which M33 lacks) is necessary to grow a massive black hole. 
A further implication is that the black hole growth in M31, M32, and the Milky Way 
was accompanied by evolution in the stellar density profiles, stellar populations, and 
dynamical structure of these nuclei such that their past appearance may have differed 
significantly from what they look like today. In short, HST observations are taking 
us to scales where understanding the central structure of galactic nuclei is intimately 
connected to the detailed interactions among their central stellar populations. 

Subject headings: galaxies: nuclei — galaxies: photometry — galaxies: structure 
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1. Introduction 

The three Local Group galaxies M31, M32, and M33 have central structures as diverse as 
their global morphologies. M31 has a double nucleus of unknown origin embedded in a bulge 
that has a relatively large core (Lauer et al. 1993). M32 has a smooth brightness profile that 
monotonically rises to extreme stellar densities but that contains no separate nuclear component 
(Lauer ct al. 1992). M33, in contrast, has an extremely compact star cluster at its center that has 
little apparent relation to the surrounding galaxy. M31 and M32 both have strong nuclear rotation 
and velocity dispersion gradients that are evidence for massive central black holes (for M31 see 
Dressier & Richstone 1988, Kormendy 1988, and Bacon et al. 1994; for M32 see Tonry 1984, Tonry 
1987, Dressier & Richstone 1988, Bender, Kormendy, & Dehnen 1996, van der Marel et al. 1997, 
and van der Marel et al. 1998), while the low rotation and velocity dispersion of the M33 nucleus 
place tight upper limits on the mass of any central dark object (Kormendy k, McClure 1993). 

One common element among these galaxies is their proximity. HST observations of all three 
systems have spatial resolution of only ~ 0.2 pc, a scale inaccessible for any other galactic system 
beyond the Local Group. Sub-parscc resolution proves crucial for even beginning to reveal the 
central properties of these galaxies. M31, M32, and M33 thus may be important "stand-ins" 
for galaxies at much larger distances. Extensive work on M31 and M32 was conducted with 
WFPC-1 and FOC prior to the spherical aberration correction. WFPC-2 observations now offer 
the opportunity to revisit these galaxies with both improved spatial resolution and dynamic range, 
with the hope of clarifying some of the issues unresolved by earlier work. 

Prior to HST, the central structure of M31 was considered to be the prototypical example of 

a dense nuclear star cluster embedded within the core of a separate bulge system. HST, however, 
showed that the M31 nucleus was double. The light that had dominated the ground-based 
images and spectra at optical and near-IR wavelengths actually came from a bright component 
(designated PI by Lauer et al. 1993) apparently physically separate and displaced by ~ 0{'5 from 
a less luminous system, P2, that itself was at the true photometric and dynamical center of 
the bulge. One obvious explanation was that PI represented the final stages of another galaxy 
being cannibalized by M31; however, this picture presented numerous difficulties related to the 
timing, dynamics, and detailed physical properties of the central components (Lauer et al. 1993). 
Tremaine (1995) suggested that the brightness peak at PI could be explained as the apogalatica 
of stars orbiting P2 in an eccentric disk. Separately, King, Stanford, k, Crane (1995) showed that 
P2 actually outshone PI at far-UV wavelengths, with perhaps some of the light associated with 
P2 having a nonthermal origin. 

The resolution of the WFPC-1 and FOC was more than adequate to reveal the important 
morphological details of the M31 nucleus. Understanding its complex form, and testing the 
Tremaine model in particular, is more likely to require higher-resolution dynamics than better 
imaging data. Nevertheless, the new WFPC-2 observations presented here have vastly improved 
dynamical range over the older data, permitting careful tests for color gradients within the nucleus 
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that might indicate interesting changes in the stellar population, as well as allowing a search for 
subtle effects of dust or other components. The data have been also taken with extremely high 
S/N and improved spatial sampling to reveal any secular morphological evolution of the nucleus 
at future epochs. 

The central structure of M32 may be typical of small elliptical galaxies in general. WFPC-1 
showed that the central stellar density of M32 may exceed lO^M© pc~^ for r < 0.2 pc,, yet its 
density profile at r ~ 5 pc, a scale accessible with HST in the Virgo cluster, is typical for a 
low-luminosity elliptical (Lauer et al. 1995). Further, M32 shows no evidence for harboring a 
separate nuclear star cluster, disk, dust, or any other structure apart from the smooth, monotonic 
rise of the overall stellar density distribution. A remaining issue is the precise form of the central 
profile at r < 0.4 pc, which is poorly constrained by the WFPC-1 images. Lauer et al. (1992) 
showed that M32 could be fitted with models that ranged from those with a small Vc = 0.37 pc 
core, to those with a 7 0.5 cusp interior to this radius, with a corresponding range in central 
density of 4 x 10^ pc"^ to 3 x lO^M© pc~^. This scale is of sharp interest as central relaxation 
and even stellar collisions are important over the age of the universe at these densities, unless the 
central dynamics are strongly dominated by a massive black hole. The new WFPC-2 observations 
of M32 have 2x the previous spatial resolution, allowing for a much improved lower limit to the 
central density of M32. Further, as for M31, the improved dynamic range allows for the detection 
of subtle color gradients that may be associated with changes in the central stellar population. 

The nuclear star cluster at the center of M33 may reach central stellar densities similar to that 
of M32 (Kormendy Sz McClure 1993). Kormendy & McClure further showed that the nucleus has 
a strong B — R color gradient in the sense that the center of the cluster is bluer than its outskirts. 
They discuss the intriguing possibility that stellar collisions at the center of the M33 nucleus have 
substantially altered its stellar population by making a young generation of blue stragglers. The 
contrast with M32 is particularly interesting. Lauer et al. (1992) argue from the lack of color 
gradients and other data that relaxation and stellar collisions have not taken place in the center 
of M32. The cross-sections for relaxation and physical collisions between stars decrease rapidly 
with increasing stellar velocity dispersion (for velocities less than the stellar escape velocities), 
due to the decreasing size of the gravitational sphere of influence of the individual stars. The 
~ 3 X 10® Mq black hole believed to be at the center of M32 (sec Kormendy &: Richstone 1995 for 
a review) would thus raise orbital speeds and greatly lengthen the relaxation and collision time 
scales in M32. Kormendy & McClure, in contrast, show that M33 has a central velocity dispersion 
of only 21 km s"-*^, thus making it highly likely that stellar relaxation and physical collisions have 
strongly modified the structure of the M33 nucleus. A basic goal of the new WFPC-2 observations 
is to obtain an improved estimate of the central density of the M33 cluster, as well as to explore 
further the centrally blue color gradient detected by Kormendy & McClure. 

We present the new images and profiles of M31, M32, and M33, followed by comparisons and 

contrasts among the three galaxies and the Milky Way nucleus in the following sections. Our 
overall conclusion is that M31, M32, and the Milky Way have had similar formation histories but 
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that M33 stands apart owing to its lack of a black hole. This difference may be related to the 
presence of a spheroidal component in the first three galaxies but the lack of one in M33. However, 
even the black-hole nuclei differ in detail, and the differences beg an explanation. An important 
part of this work is understanding the resolution limits of the new WFPC-2 data. Although the 
spherical aberration correction in the WFPC-2 optics represents a vast improvement over the 
aberrated WFPC-1 observations, the ifS'T-l- WFPC-2 PSF remains important, and its blurring 
effects must be accounted for in quantitative work by either deconvolution or modeling. 



2. Observations and Analysis 

2.1. Observational Material 

The WFPC-2 observations of M31, M32, and M33 are summarized in Table ||. The galaxies 
were centered in the high-resolution PCI CCD for all images. The exposure times were set with 
the general goal of exceeding S/N 100 in the central pixel of the galaxies, although time 
constraints prevented the F300W exposures of M31 from achieving this level. In addition, little 
if any galaxy light was seen in the F160W exposures of M31 despite the long exposure, due to 
the poor far-UV sensitivity of WFPC-2. A low gain setting (15e~/DN) was used for the F555W 
and F814W exposures to reach S/N levels high enough to measure the separation of PI and 
P2 with especially high precision. The V and I M31 exposures were also dithered by 0.5 PC 
pixels, mapping out a 2 x 2 subpixel pattern to counter the undersampling present in PC images. 
Dithered F555W images of the quadrupole gravitational lens 2237-1-0305 (discovered by Huchra et 
al. 1985) were also obtained near in time to M31, with the lens positioned at the same location 
on the PCI CCD. These observations are not presented here but can serve as an astrometric 
reference for M31 observations taken at future epochs. Deep exposures at the normal gain setting 
(7e-/DN) were also obtained in F555W and F814W for M31, and F555W for M33, to look for 
faint point sources near the centers of these galaxies. 

Basic reduction of the images was accomplished with calibration products kindly provided 
by the WFPC-2 Investigation Definition Team (IDT). The observations were obtained with 
the telescope guiding in fine-lock, and, with the exception of the dithered F555W and F814W 
exposures of M31, the multiple images in each filter sequence were generally registered to better 
than 0.05 pixels, enabling excellent cosmic ray rejection. The guiding performance for all images 
was excellent, with the jitter always being 0''004 or less. 

Cosmic ray rejection was complex for the dithered M31 images, as the 0.5 pixel shifts were 
enough to cause significant changes in the image structure that could be mistaken for cosmic ray 
events by a simple statistical examination of the pixel values. The rejection of cosmic rays was 
done iteratively; cosmic rays in a given image were identified by using an average of the remaining 
frames shifted to its position for comparison. This process was repeated a second time after initial 
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identification of cosmic rays in the comparison frames. 

The final F555W and F814W images of M31 were built from the dithered sets by simply 
interleaving the four images for each filter. The images used for analysis thus have pixels spaced 
every 0''0228, or 0.5 PC pixel. These images can be viewed as M31 convolved with the HST+PC 
PSF and a box the size of a PC pixel, but then sampled with a grid 2x finer than that of the PC. 
The F555W and F814W M31 images are thus Nyquist-sampled. 

2.2. Image Deconvolution 

Despite the dramatic improvement in image quality provided by the corrective optics installed 
in WFPC-2, PSF blurring is still important in the new data, and deconvolution or modeling of 
the images will still be required for quantitative analysis. This can be seen clearly in Figures 
H and which compare the F555W brightness profiles for M32 and M33 before and after 40 

M32 F555W Surface Brightness Profiles 

10 



11 



12 



13 



14 

.01 .1 1 

Semimajor Axis (arcsec) 

Fig. 1. — The inner 1" of the WFPC-2 deconvolved F surface brightness profile of M32 compared to 
the WFPC-2 profile measured prior to deconvolution and the deconvolved WFPC-1 profile (Lauer 
et al. 1992). Note that the deconvolved WFPC-1 profile is more accurate for r > O'.'l than the 
WFPC-2 profile prior to deconvolution. 
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M33 F555W Surface Brightness Profiles 
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Fig. 2. — The inner 1" of the WFPC-2 deconvolved K surface brightness profile of M33 compared to 
the profile measured prior to deconvolution, and the deconvolved WFPC-1 profile. Again note the 
excellent agreement of the WFPC-1 and WFPC-2 deconvolved profiles. The CHFT Bl (.B-band) 
profile of Kormendy & McClure (1993) is also shown for comparison, assuming B — V = 0.45. 



iterations of Lucy-Richardson deconvolution (Richardson 1972; Lucy 1974) to the deconvolved 
WFPC-1 profiles of the same systems. The central few pixels of the profiles as observed with 
WFPC-2 are clearly depressed compared to the same profiles after deconvolution. Indeed, in M32 
the effects of blurring persist at greater than 5% of the local surface brightness as far out as 0''5. 
In contrast, the deconvolved profiles of both WFPC-1 and WFPC-2 are in excellent agreement at 
the few percent level for r > 0''1, with the superior resolution of WFPC-2 evident in the central 
two points in the profiles. Deconvolution simulations of WFPC-1 data conducted by Lauer et al. 
(1992) had predicted that the WFPC-1 data could provide accurate brightness profiles for r > 0''1 
after deconvolution; the current deconvolved WFPC-2 profiles provide direct verification of this. 
Indeed, as can be seen in both Figures || and the WFPC-1 deconvolved profiles actually prove 
more faithful to the brightness distributions of the galaxies than do the WFPC-2 profiles prior to 
deconvolution. 

To explore the effects of deconvolution further, deconvolution simulations are presented 
in Figure ^ on the same M32 galaxy model simulated by Lauer et al. (1992) but with S/N 
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M32 WFPC-2 Cusp Deconvolution Tests 
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Fig. 3. — Surface brightness comparisons between the Lauer et al. (1992) cusp model of M32 and 
observed and deconvolved profiles. The dashed line shows the difference between the model and the 
profile as observed in a WFPC-2 image after convolution with the PSF. Note that while the largest 
differences occur for r < O'.'l, which approaches the scale of the PSF core, blurring effects are still 
visible at r ~ 0"5. The solid traces show how the profile converges to the model as a function of 
the number of Lucy-Richardson deconvolution iterations. 

appropriate to the present observations. The figure shows the residuals of profiles at various 
levels of deconvolution as compared to the input model. The dashed line shows the difference 
between the model and observation prior to deconvolution. verifying the difference between the 
WFPC-2 observed and deconvolved profiles of M32 shown in Figure |^. As can be seen, most of the 
deconvolution correction takes place in the first few iterations of Lucy-Richardson. Ten iterations 
restore the profile to about the same level of resolution as did 80 iterations on WFPC-1 images. 
After 40 iterations, the profile matches the model with errors of typically only a few percent over 
the entire profile; for r > 0^'04 (outside the central pixel), the deconvolved profiles can be taken as 
a faithful rendition of the intrinsic galaxy profile to ~ 1%. In this particular experiment, even the 
central pixel is recovered at the few percent level. As was emphasized by Lauer et al., however, 
the exact central deficit depends on the input model; the particular residual in this experiment 
should not be applied as a general correction. Again, image modeling is likely to be required at 
this stage. The final residual blurring of the central pixel is also considerably smaller in WFPC-2 



- 9 - 



than in WFPC-1, allowing the possibility of using an integrated central luminosity constraint, 
which was essentially impossible with WFPC-1. 

Evaluating the deconvolved profiles also requires understanding the effects of uncertainties in 
the PSF. We evaluate this potential source of systematic error by using different stellar images as 
representations of the PSF. Fortunately, it appears that plausible variations in the PSF produce 
both clear and benign effects in the deconvolved profile. The sole effect of plausible uncertainties 
in the core of the PSF is to modulate the brightness of the central pixel relative to its immediate 
neighbors. In particular, we find that the central pixel may vary by 0.06 mag when deconvolved 
with different PSFs. A more important problem is how the galaxy images are centered within a 
pixel, as at F555W the WFPC-2 images are undersampled. Such undersampling causes aliasing, 
which can artificially broaden sharp cusps or nuclei at the galaxy centers. We discuss the particular 
effects of this error on M32 and M33 below; but, in brief, after accounting for aliasing one can 
discriminate among various interesting possibilities for the intrinsic central brightnesses of the two 
galaxies. A final concern is allowing for the extended wings of the PSF. The wings only contribute 
to low wavenumbers in the Fourier domain; thus they mainly affect the accuracy of the overall 
slope of the profile on arcsecond scales rather than at the HST resolution limit. Nevertheless, 
because color gradients are generally very shallow in galaxies, accurate representation of the PSF 
wings in all colors is vital for measuring such gradients. 

The PSFs used for the deconvolutions were constructed from several individual star exposures 
obtained over the time span of the galaxy observations. The bulk of the images came from a 
routine monthly calibration program that imaged a standard star centered in the PC. Other 
data came from WFPC-1 IDT stellar photometry programs. For the F555W PSF, the dithered 
gravitational lens images of 2237-1-0305 provided excellent sampling of the PSF core. In all cases, 
only deep exposures that had well-exposed wings were used. The composite PSFs have total 
exposure levels of ~ 10^e~ and extend out to ~ 3"; the size of the images is 256 x 256 pixels, or 
11''6 X ll''6 (the PSFs used for the dithered M31 images have the same angular extent but double 
the sampling). 

2.3. M31 

The F555W, F814W, and F300W images of the M31 nucleus after deconvolution (with 40 
iterations) are shown in Figures ^ and |^. Note that the F555W and F814W images have been 
constructed from dithered image sets and thus have double the normal PC pixel scale, while 
the F300W image is undithered. The images show the same gross morphology that was already 
evident in the WFPC-1 images (Lauer et al. 1993); thus we will not repeat the extensive set of 
decomposition experiments that were tried by Lauer et al. (1993) and King et al. (1995). We 
do, however, present the V-band brightness profile as measured with the portion of the image 
dominated by PI excluded from the fit in Table ^ One obvious consequence of the improved 
dynamical range of the new data is the detection of strong surface brightness fluctuations in the 
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JPEG Picture 



Fig. 4. — V, I, and i7 deconvolved images of M31 (top, bottom, and top panel on next page). All 
panels are shown with an arbitrary logarithmic stretch. North is 55.7° to the left of the top of the 
pictures. The region shown is the central 9''10 x 9('10 (200 PCI pixels) centered on P2; however, 
the V and / images are 2x2 subsampled. The U image is photon-noise dominated, but the V and 
/ images are SBF dominated. 



JPEG Picture 



Fig. 5. — A "true-color" estimate of the M31 nucleus generated by using F300W for blue, F555W 
for green, and F814W for red. The intensity stretch is logarithmic. The region shown is the central 
11{'65 X 11'.'65 (256 PCI pixels) centered on P2. North is 55.7° to the left of the top of the picture. 
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y and / images, visible as strong mottling (Figure ^). Ajhar et al. (1997) have measured surface 
brightness fluctuation (SBF) amplitudes for the images presented here. In F555W, the typical 
SBF "star" has my = 25.8. At the central surface brightness in PI of 13.4, there are thus ~ 190 
SBF stars in each pixel, giving a relative noise level of 7%. In F814W, the surface brightness of 
PI is 12.0, which for SBF mi = 23.4 implies a noise level of 12%. Clearly at i/S'T resolution, we 
are running out of the very stars that trace the structures that we are trying to understand; the 
dominant error is not photon noise but the finite number of stars that make up M31 itself. SBF 
is less apparent to the eye at F300W, but this image clearly contains a number of bright stellar 
sources surrounding the nucleus. The pattern of these sources matches that seen in the King et 
al. (1992) FOC image taken at F175W; King et al. argued that these are PAGB stars. 

Another improvement provided by the WFPC-2 images is a clearer picture of the structure of 
P2 and its immediate surroundings. P2 apparently corresponds to both the center of the galaxy 
and the massive black hole inferred to exist within the nucleus (see below). The F300W image 
is particularly interesting, as the dominance of the two peaks has reversed from that seen in the 
redder colors; at 3000 A a compact blue source is strongly evident, centered on P2 (Figure 
King et al. (1995) already showed that P2 was much brighter than PI in the far-UV (at 1750 A); 
we now see that the cross-over must occur between F300W and F555W. In the WFPC-2 F555W 
image, the compact source at P2 is still evident, if greatly reduced in amplitude. Normalized 
surface brightness distributions centered on P2 are shown in Figure |^ for the three colors presented 

M31 P2 Surface Brightness Profiles 
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Fig. 6. — Surface brightness profiles of the region in M31 immediately surrounding P2. The F300W 
and F814W profiles are scaled to match the F555W profile at 0''2. 



here. The strength of the compact source in F300W is clearly evident; even at F555W the source 
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causes the profile to rise above a shallow power-law cusp at r < 0''05. Only at F814W is there no 
evidence of a source superimposed on the cusp. In passing, we note that WFPC-1 showed only 
that the F555W P2 profile consisted of a shallow cusp; the limited dynamic range of the aberrated 
observations made it exceedingly difficult to detect the faint compact central source against a 
bright background. 

The blue P2 source is compact, but not point-like — it is slightly resolved. This is evident 
in Figure |^ where we show the F555W image subtracted from the F300W image after it was 



Fig. 7. — A residual image showing the blue source at P2, creating by subtracting the deconvolved 
F555W image from the deconvolved F300W image. Note that the emission associated with P2 is 
extended. The region shown is the central d'.'lO x 9"10 (200 PCI pixels) centered on P2. North is 
55.7° to the left of the top of the picture. 

sampled to match the P2 centroid in the F300W image and scaled to match the surface brightness 
of the outer portions of the nucleus. The residual blue excess of the P2 source clearly extends over 
several pixels, while the stellar sources noted above are only 1-2 pixels in size. The F300W profile 
(Figure ^) shows that the source has a half-power radius of 0''06, or 0.2 pc. We estimate its flux 
by integrating the excess blue light within a 4 pixel radius in the difference image. The apparent 
luminosity of the source is msoo = 18.6 it 0.3. By comparing the F555W and F814W images, we 
similarly conclude that the P2 source has my = 18.7 it 0.3. King et al. (1995) raised the possibility 
that the P2 source was nonthermal emission from a low-level AGN and found it convenient to 
present their F175W luminosity in flux units. Using the zero-points given in Holtzman et al. 
(1995), we find /300 = 60 ± 20 fi3y and 7555 = 160 ± 50 ^Jy, where we have adopted Ay = 0.24 
(Burstein & Heiles 1984) and A300 = 0.42 (based on extrapolation from the 06 extinction table 
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of Holtzman et al.). These fluxes can be compared to /175 = 7 iJ,3y observed by King et al. The 
physical nature of the P2 source is discussed in § ^.4 

The new data confirm the conclusion reached by Lauer et al. (1993) that P2 corresponds 
to the photocenter of M31. The location of isophote photocenters as a function of V surface 
brightness is shown in Figure |8|. Centers were calculated by measuring the "center of mass" of all 

M31 F555W Nucleus Isophote Centers 
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Fig. 8. — The RA and DEC offsets of the isophote centers with respect to P2 are shown as a 
function of F555W surface brightness. The first point plotted near PI is at nv = 13-1; succeeding 
points step 0.1 mag in surface brightness. The solid points mark isophotes of = 13.5, 14.0, 14.5, 
15.0, and 15.5. All isophotes with fiy ^ 14.5 are centered within 0''1 of P2. 



pixels within ±0.1 mag of a given isophote level. Isophotes with 15 < fiy ^ 16 all have centers 
within O^'l of P2. Indeed, the centroid of pixels with /uy = 15.5 it 0.3 is only ~ 0''05 from P2. As 
can be seen in Figure which shows the central brightness profile of M31 after excluding the area 
strongly affected by PI, isophotes of this brightness are ~ 2" in radius and thus define the outer 
portions of the nucleus (that is the portion of the brightness profile that rises above the bulge 
profile). Isophotes at larger radii are affected by dust but also appear to have centers within 0"1 of 
P2. Notably, the offsets seen in these cases are perpendicular to the P1-P2 vector (see the lowest 
surface brightness points in Figure ^); their projection on the P1-P2 vector is thus still close to P2. 

Color variations appear to show that the center of the nucleus (the region within 1" of P2) 
has a slightly different stellar population than the outer portions of the nucleus and surrounding 



bulge. Figure 10 shows the ratio of the F555W and F814W images. PI is clearly redder in 
V — I than the surrounding galaxy, but the region of red excess also extends into the region 
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surrounding P2 on the opposite side from PI. The central 0''4 x 0''4 region of PI has V — I = 1.41, 
while the color in a similar-sized aperture displaced 0''4 from P2 in the "anti-Pl" direction has 
V — I = 1.38. The bulge and outer portion of the nucleus are bluer still with V — I = 1.34 for 
l" < r < 5" . On the other hand, PI and its anti-region are bluer than the bulge in F300 — V, 
with PI having F300 - V = 2.29, the anti-Pl part of P2 having F300 - V = 2.18, and the 
isophotes with 1" < r < 5" having F300 — V = 2.39. In passing, we note that this result is in 
apparent contradiction to the results of King et al. (1995), who found rec/tfer-than-bulge colors 
in F175W - F555W from pre-COSTAR FOG and WFPC-1 data in the region outside r ~ 1". 
For comparison, P2 itself has V — I = 1.28, and F300 — V = 1.29. For convenience, these color 
measurements are summarized in Table |2|. 

Semimajor Axis (parsec) 




.01 .1 1 10 

Semimajor Axis (arcsec) 

Fig. 9. — The F555W brightness profile of the M31 nucleus. The portions of the image dominated 
by PI have been excluded from the profile measurement. 



That the center of the nucleus is redder than the bulge inV — I yet bluer in F300 — V may be 
reasonable; such combinations can in principle be created by simultaneous changes to both stellar 
population age and metallicity. However, an indisputable conclusion from the color difference is 
that the central nucleus contains stars with an origin different from the surrounding bulge. This 
may support the Tremaine (1995) eccentric disk model of the doubly-peaked nucleus; the disk 
would consist of material that has fallen into the nucleus after the central structure of M31 was 
already organized. Additional support for the Tremaine disk may come from the morphology 
of the red-excess region evident in Figure 10, given that it clearly extends to both sides of P2; 
PI would simply be a local density enhancement in a population with a more extended spatial 
distribution. The intermediate color of the anti-Pl part of P2, that is a color between that of 



- 16 - 



m 

CD 

■v 

o' 

c 
—1 

CD 



Fig. 10. — The ratio of the M31 F555W and F814W deconvolved images. The scale is magnitudes, 
with white corresponding to V — I = 1.1, and black to V — I = 1.6. The region shown is the central 
d'.'lO X g^'lO (200 PCI pixels) centered on P2. North is 55.7° to the left of the top of the picture. 
P2 is evident as the white peak at the very center of the image, surrounded by a region of redder 
than average color. The noisy appearance of the image is due to the incipient resolution of M31 
into stars. Tick marks point to the centroids of PI and P2 

PI and the bulge, may reflect underlying dilution of the disk by the bulge. It is also noteworthy 
that the intensity minimum between PI and P2 that defines the double-peaked appearance of the 
nucleus has V — I = 1.34 it 0.02, a color bluer than that of PI but which matches the bulge. 



2.4. M32 

The F814W image of the M32 nucleus after deconvolution (with 40 iterations) is shown in 



Figure |lT| a color image is shown in Figure 12. Surface photometry profiles are presented in 
Figure |^ and Table ^. The surface photometry was measured with the algorithm of Kent (1983), as 
modified by Lauer (1985) for high-resolution applications. The algorithm uses sinc-interpolation 
and Fourier harmonics to measure the brightness distribution along isophotes at every integer-pixel 
step in semimajor axis length. 

The most impressive feature of M32 is its incipient resolution into stars, as was also seen 
in M31. The surface brightness fluctuations in F814W are especially strong; as in M31, they 
dominate the profile error budget. Ajhar et al. (1997) have measured SBF amplitudes for the 
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Fig. 11.— The F814W deconvolved image of M32. North is 120.0° to the right of the top of the 
picture. The region shown is the central 23''3 x 23''3 (512 PCI pixels). The stretch is logarithmic. 



JPEG Picture 



Fig. 12. — A "true-color" image of M32 formed from the F555W and F814W deconvolved images. 
North is 120.0° to the right of the top of the picture. The region shown is the central 23''3 x 23''3 
(512 PCI pixels). The stretch is logarithmic. 
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images presented here. In F555W the typical SBF "star" has my = 25.4. With a central surface 
brightness of 10.96, or 17.67 per PC pixel, there are ~ 1, 200 SBF stars in the central pixel, giving 
a relative noise level of 2.8%. In F814W, //q = 9.79 and m/ = 23.0, implying a noise level of 5.0%. 
These errors are close to the deconvolution uncertainties in the central pixel. 

Aside from the strong SBF pattern, we find no evidence for dust, disks, or any other structures 



present in M32 down to levels of a few percent in local surface brightness. Figure 13 shows a color 



ratio image made by dividing the V image by the / image, and the residual map resulting from 






Fig. 13. — The left panel shows the color-ratio image of M32 formed by dividing the deconvolved 
F814W image into the F555W image. The right panel show a residual map formed by subtracting 
a model reconstructed from the M32 F555W surface photometry profiles from the original image. 
North is 120.0° to the right of top in both panels. The region shown is the central 9"10 x 9"10 (200 
PCI pixels). The stretch is ±1 mag is surface brightness about the mean residual. 

subtraction of models reconstructed from the F555W profile. All such comparisons are flat on 
average (although the color-ratio image shows dark "pits" due to the red color of barely resolved 
bright stars). The isophote shape profiles additionally show no evidence for any structures not 
accounted for by elliptical isophotes. Figure |l^ shows the A4 isophote shape parameters for both 
F555W and F814W, which give the relative fourth harmonic deviations of the isophotes from 
perfect ellipses. Both profiles are noisy due to strong surface brightness fluctuations but show no 
evidence for the isophotes to be significantly "boxy" or "disky" (which would give significantly 
negative or positive A4, respectively). This verifies the conclusions of Lugger et al. (1992) and 
Lauer et al. (1992) and stands in contrast with the claim by Michard &: Nieto (1991) that the 
inner isophotes of M32 were disky, having positive ^4. Michard & Nieto specifically argued that a 
disk was visible in their reddest images and thus might be comprised of AGB stars. The F814W 
profile, while noisier than F555W, shows no greater evidence for a disk. It is also noteworthy 
that the ellipticity of M32 remains constant throughout the inner portion of the galaxy, having 
€ = l-b/a = 0.274 ± 0.003 for r > 0''1 (Figure |l|). Lauer et al. (1992) argued that the isophote 
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Fig. 14. — Isophote shape profile of M32 for both the F555W and F814W deconvolved images. 
Positive A4 corresponds to "disky" isophotes; negative A4 corresponds to "boxy" isophotes. 
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Fig. 15. — Isophote ellipticity profile of the deconvolved F555W image of M32. The dashed line 
indicates the average ellipticity for 0('l < r < 10". 
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ellipticity in M32 stayed constant into the resolution limit, based on deconvolution simulations 
with WFPC-1; the constant ellipticity is shown directly here with WFPC-2 data. 



As noted in §2.2, the new WFPC-2 observations of M32 are in excellent agreement with the 
WFPC-1 F555W profile of Lauer et al. (1992) for r > O'fl but rise above the old WFPC-1 profile 
at smaller radii due to the improved resolution of WFPC-2. The new data allow us to resolve 
some of the ambiguity in the central structure of M32 due to uncertainties in the WFPC-1 profile. 
Lauer et al. (1992) showed that the central structure of M32 was well bracketed by a "core" model 
of the form: 

^(0=^0 l + 7v7T7^ , (1) 



where in the F-band Iq corresponds to /ig 



*v 0"166y 
10.99, and a "cusp" model of the form: 



I(r) = In 



-0.53 



1 + 



0^'378 / 



-0.375 



(2) 



where Iq corresponds to fio 
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13.04. Both of these models are shown in Figure IC. All but the 
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Fig. 16. — The M32 V surface brightness profile and models. The "cusp" and "core" models were 
introduced by Lauer et al. (1992) and are discussed in the text; the Nuker model is new to this 
work. The horizontal lines bracketing the central point of the deconvolved profile show the central 
brightnesses of the core and cusp models accounting for aliasing; the Nuker model corrected for 
aliasing is not shown, but fall slightly below the cusp model. The Nuker model now gives the 
best fit. The inset figure shows the mass-density profile implied by the Nuker model, assuming 
M/Ly = 2.0 in solar units. 



central point of the V profile fall on the cusp model. The core model fits less well, rising above 
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the observed profile for O^'OS < r < 0''3 and falling below the central two points. Since the 
deconvolution simulations show that the profile is accurate to ~ 1% at all but the central point, 
the cusp model can be taken as an excellent description of the M32 profile for r > 0''04, or 0.3 pc. 

Understanding the light distribution of M32 at smaller radii clearly depends on how well we 
can use the central point in the brightness profile as a constraint. If we take the effective radius of 
the central pixel as 2~^^'^ of the pixel width, or 0''016, then it is offering information on the central 
0.1 pc of M32. Superficially, the central point of the F profile falls below the cusp model by 0.34 
mag. It appears that this deficit is not due to residual blurring from incomplete deconvolution. 
Simulated observations of the cusp model presented in § |2.2| show that even the central point in 
the deconvolved F555W profile can be recovered to ~ 1% with an accurate PSF. Further, the 
deconvolved image matches the observed image extremely well when it is convolved back to the 
resolution of the latter; if the deconvolution were incomplete, its central brightness would be 
too low to match the observation after this operation. Fitting the observed image with the PSF 
convolution of the deconvolved image plus a scaled PSF to make up for any central shortfall in the 
deconvolved image shows that the central pixel is likely to be underestimated by only 0.01 it 0.01 
mag. Systematic errors in the PSF are important and will cause corresponding profile errors larger 
than this amount. Using a different PSF that was less well sampled than the dithered F555W 
PSF used for the present results gives a deconvolved profile 0.06 mag brighter in the center. As 
noted above, the intrinsic error in the central brightness itself is 0.03 mag due to the finite number 
of central stars, in good agreement with the 0.04 mag of noise measured in the central 5x5 pixel 
patch by comparing the image with a model reconstructed from the photometry. Adding the PSF 
error and surface brightness error in quadrature gives the total error in the central pixel of 0.07 
mag. 

The most serious difficulty with measuring the central brightness is that at F555W the PC 
images are undersampled by ~ 2x, and are thus aliased. Aliasing causes spatial information 
above the Nyquist sampling frequency to contribute erroneously to structure at frequencies below 
the sampling frequency. In simpler terms, it causes irreversible broadening of features. This in 
fact appears to be happening at the very center of M32. Unfortunately, the photocenter of M32 
appears to be significantly displaced from a pixel center. The CCD-row offset of the photocenter 
is only 0.15 pixels, but the CCD-column offset is 0.46, causing significant broadening of the M32 
nucleus in this direction. To understand the quantitative effect of aliasing, we constructed an 
image of the M32 cusp model faithful to the orientation and centroid of the observation. The 
model image was constructed with 3x subsampling of the PC pixels, with additional subsampling 
of the central portion of the cusp itself. When a profile was measured from the synthetic image 
positioned and binned to match the PC sampling, the value of the central pixel was fainter by 
0.21 mag than the input model, leaving the the central difference between the cusp model and 
observed profile at only 0.13 it 0.07 mag, no longer highly significant. A similar treatment of the 
core model shows that, once aliasing is accounted for (and the core model is rescaled in brightness 
to match the observed profile at ~ 0^'2), the central observed point exceeds the central brightness 
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of the core model by 0.21 it 0.07 mag. The ahased central points of both models are shown in 
Figure |l^ as horizontal lines bounding the observed central point. 

The conclusion is that the center of M32 is significantly more concentrated than the core 
model. Despite aliasing corrections, the central point still falls below the cusp model, although 
with low significance. Lauer et al. (1992) concluded that the asymptotic slope of the cusp model 
has 7 = 0.53 it 0.05. Taking the corrected central point as accurate implies 7 = 0.43 it 0.07 for 
< 0''04. A Nuker-law fit of the form (Lauer et al. 1995) 



I(r) = 2('3-^)/"Ife ( ^ 



7 r 



, r 
1+1 - 

n 



(7-/3)/a 



(3) 



to the y-profile for r < l" gives a = 1.39 ±0.82, (3 = 1.47 ±0.16, 7 = 0.46 ±0.14, h corresponding 
to fib = 12.91 ± 0.31, and = 0^'47 ± 0''15; this fit is shown in Figure 16; note that this model 



is essentially identical to the cusp model at all but the central pixel, where it predicts slightly 
less luminosity. In passing, we note that the large errors in the fit parameters include their high 
covariance; the overall form of the model is tightly constrained over the domain fitted (Byun et al. 
1996). Taking the Nuker model as the best representation of the profile implies that the central 
mass density in stars (Figure 17)is at least 10'' Mq pc~^ at the HST resolution limit (assuming 
M/Ly = 2.0 in solar units. Ay = 0.24, and 770 kpc distance). Clearly M32 remains "cuspy" at 
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Fig. 17. — Stellar mass density profiles of M31, M32, and M33 based on Abel inversion of the 
F555W brightness profiles, assuming stellar population M/Ly = 6.5, 2.0, and 0.4 in solar units for 
the three galaxies, respectively. For M31 and M33 direct inversions of the deconvolved brightness 
profiles are shown; for M32, inversion of the best-fitting Nuker model is shown. 



the resolution limit of HST. Answering the question of whether we are seeing the cusp start to 
flatten, however, should be addressed by obtaining new images with sub-pixel dithering. 
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At larger radii, the WFPC-2 observations place excellent constraints on color gradients that 



were not possible with WFPC-1. As noted in §2.2, even with WFPC-2, the effects of the PSF still 



extend to surprisingly large radii, and accurate color gradients require deconvolution and a PSF 
accurate on large spatial scales. In the case of M32, it appears that the lack of color gradients at 
radii observable from the ground (for U — R, Michard & Nieto 1991; and B — R, Lugger et al. 
1992) continues into the center at least in V 
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Fig. 18.— V - 
O'fl <r < 10". 



/ and U — V profiles of M32. The dashed lines show the mean color gradients for 



to theV - I trace for O'.'l <r < 10" gives V - I = (1.237 ± 0.002) - (0.009 ± 0.005) log(r/l"), 
which is essentially flat. The color for r < 0"1 becomes tantalizingly bluer. However, the intrinsic 
color error of the central point is 0.06 mag from SBF alone, and 0.10 mag when PSF errors are 
included — the apparently bluer color is not significant. As a check on this result, we obtained 
F336W PC exposures of M32 from the STScI archive and processed them as we did the present 
observations. The resulting U — V color gradient is also shown in Figure |l^; the fitted gradient is 
U -V = (1.216 ± 0.004) - (0.023 ± 0.008) log(r/l"). Further, Cole et al. (1998) found no far-UV 
(1600 A) versus V gradient for r < 0"8. The U — V colors should be a more sensitive test for hot 
stars at the center of M32 than V — I; as the U — V color becomes redder towards the center, we 
find no evidence that M32 harbors any population similar to that at the center of P2. From the 
errors, we estimate that any blueing trend in the central pixel of M32 is no more than 15% of the 
U — V blueing seen within r = 0.2 pc in M31. 
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2.5. M33 

The F555W image of the M33 nucleus after deconvolution is shown in Figure |l^; a color image 
is shown in Figure Surface photometry profiles are presented in Figure ^ and Table |5[ Both 




Fig. 19. — The F555W deconvolved image of M33. North is 113.4° to the right of the top of the 
image. The region shown is the central 4''55 x 4''55 (100 PCI pixels). The stretch is logarithmic 
and covers 10 magnitudes in surface brightness. 

the images and surface photometry plots show vividly that M33's nucleus is extremely compact. 
While the central F555W surface brightness is brighter than 11 mag/D arcsec, the brightness 
profile falls rapidly as a power law with slope 7 = 1.49 it 0.09 for the first semi-decade in radius 
(O^'OS < r < 0"25), steepening to 7 = 1.90 ± 0.03 over the next semi-decade (0"25 < r < 0^8). At 
r ~ 2", the nucleus profile begins to fall below the disk, having decreased to > 18.5, or ~ 10^ 
dimmer than the central value. At r ~ 1", the nucleus may be showing signs of incipient resolution 
into stars, but as the disk begins to dominate, it is difficult to associate any of the point sources 
at larger radii with the nucleus itself. 

We first consider the WFPC-2 constraint on the central luminosity density. Kormendy 
&; McClure (1993) argued that the M33 nuclear core radius must have < O'-'IO, based on 
seeing-convolved models. As shown in Figure the present observations match the best Kormendy 
& McClure i?-band profile for r > 0''3 (assuming B — V = 0.45), but climb to surface brightnesses 
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Fig. 20. — A "true-color" image of M33 formed from the F555W and F814W deconvolved images. 
North is 113.4° to the right of the top of the image. The region shown is the central 11''7 x 11''7 
(256 PCI pixels). The stretch is logarithmic. 

over an order of magnitude brighter at small radii. Near the HST resolution limit, the M33 
profile is still climbing steeply. If we take the deconvolved profile as accurate excluding the central 

point, then the first and second points out give 7 = 1.44, consistent with the slope measured at 
O^'OS < r < 0'.'25. The central point a priori is most likely to be a lower-limit on the true central 
surface brightness. Dcconvolution experiments on simulated M33 profiles show that we can recover 
the central intensity at the few percent level, provided that the PSF is accurate; however, as 
with M32, aliasing of the images is important. The photocenter of M33 in the F555W image is 
displaced 0.41 pixels in the CCD row direction and 0.30 pixels in the column direction from the 
center of a pixel. 

As with M32, we constructed models with steep central power-law profiles on a subpixel grid, 
and then binned the models to match the observations. For cusps with 7 ~ 1, the aliasing of the 
M33 image depresses the central intensity from its true value by over 0.2 mag. The procedure for 
estimating the central brightness of M33 was to assume that its intrinsic profile continues inward 
as a constant power-law for r < 0"04 and construct models for a range of slopes that produced 
central decrements bracketing the observations. The best fitting model gives 7 = 0.82 it 0.17 after 
accounting for aliasing; the uncertainty reflects an adopted error in the observed central intensity 
of 0.07 mag due to PSF uncertainties and surface brightness fluctuations. This result suggests that 
the M33 nucleus profile is still becoming shallower near the HST resolution limit, continuing the 
slow decrease in logarithmic slope seen over 0''08 < r < 0''8. Inward extrapolation of the 7 = 1.49 
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slope seen just outside the center predicts an intensity decrement between the central and next 
pixel in the profile of 0.94 magnitudes, compared to the observed decrement of only 0.59 mag. 

In contrast, Kormendy & McClure (1993) fitted their profiles to a form that included an 
"analytic" core, that is, one that requires 7^0 smoothly as r ^ 0. In this spirit, we also modeled 
the central portion of the profile with the form: 

Fits to the 1^-band image give a = 0.6 it 0.1 PC pixels. Defining rc as the half-power point implies 
Tc = 1.24a, thus Tc = 0''034 ± 0''006, or 0.13 ± 0.02 pc for a distance to M33 of 785 kpc. This form 
is no more favored than the constant 7 = 0.82 cusp; thus this core size should be regarded only 
as an upper limit. We note that Massey et al. (1996) claimed r^. = O^'OS from the present images 
obtained from the STScI archive; however, they neglected to account for the effects of the PSF or 
aliasing. 
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Fig. 21. — Central color gradients in M33. The upper panel shows the difference between the V 
and / profiles as a function of radius (uncorrected for reddening). The bottom panel compares the 
present V — I profile degraded to CFHT seeing (solid line) to the B — R profile of Kormendy &: 
McClure (1993) (open circles). 
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Abel inversion of the M33 brightness profiles implies that its nucleus reaches central mass 
densities of at least 2 x 10^ pc~^ (Figure 17). While the peak luminosity density of M33 
exceeds that of M32, Kormendy & McClure (1993) conclude that M33 has central M/Ly < 0.4, 
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a value ~ 5x less than that of the stellar population of M32. M32 thus reaches much higher 
mass-densities on the scales probed by HST. 

The structural analysis was performed on the y-band image, as it has higher spatial resolution 
than the /-band image. Comparison of the F and / profiles reveals a strong central color gradient, 
which appears to be consistent with the B — R color gradient observed by Kormendy & McClure 
(1993) at CFHT resolution. Figure 21 shows the V — I color profile obtained from the two 
deconvolved profiles. The color appears to be constant at V — I = 1.00 it 0.01 for r > O'.'S but 
becomes bluer, with V — I = 0.85 it 0.02, closer to the center. The low surface brightness of the 
nucleus at r ~ l" yields large errors with the fine-resolution brightness profiles used to generate 
the color profiles; using larger bins reduces the noise but still shows that V — I is essentially 
constant for r > O'.'S. In passing, we note that even — / = 1.0 is significantly bluer than the 
typical y — / ~ 1.3 for the central regions of elliptical galaxies (Lauer et al. 1998), consistent with 
the well known blue color of the M33 nucleus (van den Bergh 1991). These colors are as observed, 
without reddening corrections. For Ab = 0.18 (Burstein & Heiles 1984), A{V — I) = 0.05, giving 
the central intrinsic V — I = 0.80 and V — I = 0.95 at larger radii. 

Kormendy & McClure find A{B — /?) ~ 0.4 over 0''2 < r < l^'O. This gradient is larger in 
amplitude and occurs at larger radii than the V — I gradient, but their data were not corrected 
for seeing. When the present data are degraded to the best CFHT seeing of 0^'44, the new V — I 
profile appears constant for r < 0^'2 but shows an essentially linear gradient for 0^'2 < r < 1''0, in 
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Fig. 22. — Isophote ellipticity profile of the deconvolved F555W image of M33. 



good agreement with the general shape of the Kormendy & McClure B — R profile (Figure 21). 
The present results thus confirm the Kormendy & McClure result that the nucleus becomes bluer 
towards the center. The larger A(i? — R) compared to IS.{V — I) is consistent with the central 
population containing hotter stars, as discussed in the next section. 

The surface photometry of M33 also includes information on the shape of the nucleus. 
Inspection of the images shows that the nucleus appears to be slightly flattened, which is confirmed 
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by the isophote ellipticities plotted in Figure 22. The average isophote shape over the central 
arcsecond has e = 0.16 it 0.01, with the nucleus becoming somewhat rounder towards its center. 
The position angles of the isophotes also appear to be constant, with average PA = 17° it 2°. 
As discussed later, the relaxation time in the M33 nucleus is extremely short; the ellipticity of 
the nucleus is probably more likely due to rotation rather than an anisotropic velocity ellipsoid. 
Kormendy & McClure (1993) find a small velocity gradient of 6 it 2 km s~^ arcsec"^ (with slit 
PA = 23°), suggestive of rotation, which can be compared to the central velocity dispersion of 
21 lb 3 km s~^. The slightly flattened shape of the nucleus appears to be consistent with this 
modest amount of rotation. 



3. Discussion 
3.1. M31 

A key advance provided by the WFPC-2 images of M31 is the detailed look at the color and 
morphology of P2. As noted above, King et al. (1995) first observed the far-UV dominance of P2; 
the present data show that P2 dominates at 3000 A as well and is much bluer in V — I than PI, 
the rest of nucleus, and the surrounding bulge. As noted by Lauer et al. (1993) and confirmed 
here, P2 closely corresponds to both the dynamical and photocenter of M31, and thus is likely 
to coincide with the 3 x 10'' black hole inferred to exist at the center of M31 (Kormendy &: 
Richstone 1995). King et al. were thus tempted to consider that the UV flux of P2 was optical 
emission associated with a low-level AGN. In support of this hypothesis they cited the radio 
source observed by Crane, Dickel, & Cowan (1992) and the X-ray source observed by Trinchieri & 
Fabbiano (1991), both coincident with the center and both having flux densities that bracketed the 
UV luminosity of P2. At the same time. King et al. realized that this picture probably required 
P2 to be a point source in the UV, an issue that was difficult to test with the aberrated FOC data. 

The present images show that P2 actually has a half-light radius of ~ 0.2 pc at 3000 A and 
is clearly extended in V as well. This result suggests the more prosaic possibility raised by King 
et al. that there is tightly bound cluster of early-type stars associated with P2. The range of 
fluxes from 1750 A to ^ can potentially be explained by a population with Tg ~ 10^ K, which 
corresponds to late B to early A-type spectral class. This is in good agreement with the FOC 
F175W and F275W flux measurements by Brown et al. (1998), who estimate Tg = 1.2 x 10"^ K, for 
P2. The smooth appearance of the P2 cusp implies that the stars are individually faint compared 
to the total luminosity; a turnoff at My 0, corresponding to ~ 10"^ K, would be suitably 
fainter than the total P2 blue source magnitude of My ~ —5.7. 

Although a cluster of warm stars might explain the blue center of P2, the origin of such a 
cluster is somewhat problematic. One obvious possibility is that the cluster consists of recently 
formed stars. The presence of gas and dust in the bulge of M31 would provide raw material for 
such stars. However, at 0.2 pc, the nominal core of the blue population, the local circular velocity 
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about the 3 x IO^Mq black hole is ~ 800 km s~^, and the rotation curve is Keplerian with strong 
shear; it is difficult to picture even a disk of dense gas here as a favorable site for recent star 
formation (we discuss this issue further in §^. 

It may be possible that stellar collisions play a role in explaining the blue population. King et 
al. (1995) briefly considered this picture but dismissed it without quantitative discussion; however, 
the new stellar density profile suggests that this idea merits closer attention. In general, the mean 
time for any star to suffer a collision in a Maxwellian distribution of identical stars is 

16V^noT^(l + G)j~\ (5) 

where n is the stellar number density, a is the 1-d stellar velocity dispersion, m^, is the stellar mass, 
r^, is the stellar radius, and = Gm^,/(2(T^r*) = (V^^'^ /2a)'^ is the "Safronov number" (Binney 
& Tremaine 1987). Abel inversion of the F555W profile centered on P2 implies that, near the 
i/^T resolution limit, the P2 stellar mass density has the form p ~ 1.3 x 10^(0.1 pc/r)^MQ pc~^, 
given the bulge M/Ly = 6.5 (Kormendy 1988), Ay = 0.24, and 770 kpc distance (Figure |l^. For 



a = 7.2 X 10^(0.1 pc/r)^/^ km s~^ implied by the nominal black hole the mean collision time 
for solar-type stars is tc ~ 6 x 10^''(r/0.1 pc)^/^yr (this calculation is summarized in Table P). 
For a stellar population with lifetime lO^'' yr, roughly 16% of main sequence turnoff stars would 
collide, and thus roughly 2x 16% of the total stellar luminosity (see §^.2|) would be radiated from 
blue-straggler-type stars of higher surface temperature. Within the accuracy of this calculation, 
this is not far from the 20% excess light contributed by the blue population of P2 in the V-band. 

On the other hand, the stellar density profile presented in Figure |l^ presumes constant 
M/Ly, which will not to be valid if P2 has been modified by stellar collisions; the central density 
will be overestimated and thus the collision time underestimated. The central blue part of P2 is 
clearly visible in the density profile as a distinct inflection; naive inward extrapolation of the outer 
part of the P2 cusp yields a more modest p ~ 4 x IO^Mq pc~^ and tc ~ 2 x lO^^yr at 0.1 pc. With 
this collision time, the contribution of blue light would be only ~ 10%. Furthermore, the 10^ K 
effective temperature of the warm component, corresponding to a mean spectral type earlier than 
AO, is hotter than blue stragglers seen in globular clusters, which are typically late A or F stars. 
It is not clear that single collisions of two main-sequence turnoff stars can yield stars this hot 
(Bailyn 1995). Hotter, more massive merger products could of course be produced by collisional 
runaway caused by the sinking of massive merged stars to the center. However, the small fraction 
of light involved in the P2 blue excess suggests that no such runaway has occurred. Moreover, the 
dynamical relaxation time at 0.1 pc in M31 is ~ 10^^ yr on account of the high velocity dispersion 
near the black hole, further evidence against such a runaway. Finally, there is no guarantee that 
massive blue stragglers would actually form from the energetic stellar collisions in M31, as the 
star-star impact velocity at 0.1 pc is ~1800 km s~^ (assuming a 1-d velocity dispersion of 700 km 



^For an r "^''^ density cusp surrounding a black hole, a — {V^ /?>) — \J2GM, / (5r), where M, is the black hole 
mass (Young 1980). A cusp of this type is not strictly correct for P2, but it is an excellent match to M32. 
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s ^ and an isotropic velocity ellipsoid). This exceeds 2x the main-sequence stellar escape velocity 
of ~ 600 km s~^, the point where major mass loss sets in in stellar collision models (Benz &: Hills 
1988). 

Note, however, that all preceding arguments implicitly assume a single hot stellar population 
with a more-or-less isotropic velocity ellipsoid. In §^ we revisit these estimates under the 
assumption that a cold disk of stars is present in the nucleus and find that the collision time in 
such a system could be shorter and more favorable. Such a cold disk is a key assumption in the 
Tremaine (1995) eccentric-disk model of the nucleus. Evidence in favor of it comes from both 
the bluer population at the center of P2 and the red color of the inner portions of the nucleus. 
The former evidence favors young stars and suggests that significant infall of gas to the center 
of M31 may have occurred well after its initial formation. The latter evidence indicates that the 
population of the nucleus differs significantly from the bulge, again supporting its formation as 
a separate event, the color map is much more symmetrical than the intensity distribution. The 
color uniformity makes it more difficult to accept PI as a separate stellar system grafted onto the 
nucleus rather than as a simple rearrangement of the basic nuclear stellar population. 



3.2. M32 

The previous WFPC-1 observations of M32 admitted a wide variation in possible central 
structure for r < 0"1 (Lauer et al. 1992). The new observations show that the surface brightness 
of M32 is still rising as a 7 ~ 0.5 cusp interior to r < 0.1 pc, reaching p > 10^ Mq pc~^ in stars 
alone. The corresponding lower limit from WFPC-1 was only ~ 4 X 10^ pc'^. At the same 
time, the WFPC-2 images extend to greater sensitivity and smaller radii the conclusion that the 
inner structure of M32 is featureless aside from the inward continuation of a strong cusp. The 
center of M32 has no detectable disk, dust, boxy isophotes, or isophote shape changes. The V — I 
and U — V color gradients are nearly flat, showing no dramatic change in the stellar population as 
one approaches the postulated central massive black hole. 

The recent FOS observations of van der Marel et al. (1997) show that the velocity dispersion 
continues to rise into the center at HST resolution, thus the rapid relaxation implied by the Lauer 
et al. (1992) "core model" without a black hole can be ruled out. In general, the relaxation time 
scale is 

_ 0.34(j3 

~ G2m,pln(0.4iV)' ^ ^ 

where p is the stellar mass density and A*" is the number of stars in the system (Binney & 
Tremaine 1987). The present observations show that the central mass density profile has the form 
p{r) « 7 X 10^ (r/0.1 pc)"^/^M0 pc'^. For a{r) = 242(r/0.1 pc)"^/^ km s"^ as implied by the 
3.4 X 10^ M0 black hole mass of van der Marel et al. (1998), t,. ~ 3 x 10^ yr throughout the center 
of M32 (assuming N ~ 10^), too long to have a significant structural effect. 
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Lauer et al. (1992) argued that a black hole of this magnitude would greatly increase the 
collision time scale over that calculated for M32 without a black hole; however, they also noted 
that significant collision effects might still be expected at scales poorly probed by the WFPC-1 
observations. From equation (^) and the a and p profiles presented in the previous paragraph, we 
estimate tc ~ 2 x lO^'^ yr at O.lpc, increasing outwards roughly linearly with radius (see Table 

This new central collision time scale in M32 may be low enough to imply detectable changes 
in the spectral energy distribution if the stellar population is old. We assume that collisions 
acting on a base population of turn-off stars of age tto and mass Mto are creating blue stragglers 
of lifetime ts ~ 10^ yr and mass Ms- Let the bolometric luminosity of the blue stragglers be L5 
and that of the turn-off stars be Lto- If the system is older than the blue-straggler lifetime, then 
creation of blue-stragglers from collisions, rito/tc, is in rough equilibrium with their death rate, 
ns/ts- Assume further that the bolometric luminosity of a star is simply L ~ [M^c^ /t^). Then 
ns = ntoits/tc)i and the ratio of bolometric luminosities emitted in steady state would be 

L5 ^ fMs_\ ftjo\ 2 X ( \ 
Lto \MtoJ\tJ UxlQiOyr/ 

If the M32 central stellar population is old, i.e., tto ~ lO^*' yr, the flux ratio is of order unity, and 
the predicted fraction of emitted blue flux would exceed that from P2 in M31 by a 5x (cf. §|3.1|. 



In fact, the color gradient discussion in §2.4 failed to find any blueing in the central pixels of M32, 



limiting it to at most 15% of that seen in M31. 

Some easing of this dilemma would occur if the stellar population were young. Estimates 
based on H/? put the light-weighted central stellar age near 2 x 10^ yr (Gonzalez 1993, Trager et 
al. 1998), which would reduce the predicted blueing effect by a 5x. However, the predicted blueing 
would still be comparable to that seen in M31, as its estimated age is also rather young (5 x 10^ 
yr, Gonzalez (unpublished)) and the ratio Uo/ts is similar in the two galaxies. Furthermore, the 
bulk of the population could be older than the lightweighted age, and we have also neglected 
the extra collisions that would occur during the post-main sequence phases. On the other hand, 
collision speeds may be high enough to prevent efficient formation of blue stragglers (see § (3.1] ). In 
sum, the stubborn refusal of M32 to reveal any departures from utter normalcy at small radii may 
be stressing simple theories of its central structure, a problem that deserves closer scrutiny with 
realistic multi-component evolving dynamical models. 



3.3. M33 

The key to the M33 nucleus may be that, unlike M32, it lacks a massive central black hole 
to stabilize its structure. Kormendy &: McClure (1993) have emphasized their impressively low 
upper limit of only M, < 5 x 10^ Mq for a central black hole in M33; with our smaller limit on the 
core radius, we can reduce this further to M, < 2 x W^Mq. With a central velocity dispersion of 
only o" = 21 km s~^ but a central density nearly as high as that in M32, the stellar collisions and 
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relaxation, in particular, will be much more important for M33. Stellar collisions indeed may be 
the explanation for both the relatively blue color of the M33 nucleus and its apparent population 
of young stars (van den Bergh 1991). 

Hernquist, Hut, &; Kormendy (1991) emphasized the short central relaxation time in the 
M33 nucleus and the likelihood that it has undergone core collapse, based on a preliminary 
analysis of the Kormendy & McClure (1993) observations. The present observations imply an 
even shorter relaxation time than estimated in either work; from equation (^) and assuming 
a = r^^^T^GpQl^ = 21 km s^^ for M33, 

The implied core-collapse time is only tec ~ 1 x 10^ yr, taking t^c ~ 330tr (Cohn 1980); it is thus 
highly likely that the nucleus has already undergone core collapse. 

Hernquist et al. (1991) were pessimistic that the blue color of the nucleus could be explained 
via blue stragglers generated in stellar mergers but did not present a quantitative argument. 
Kormendy & McClure (1993), in contrast, were intrigued that stellar collisions might indeed be 
important. Simple calculations of the central stellar collision time scale suggest that collisions are 
highly likely over the age of the nucleus. For the case with o" << GMjr^^ equation (^) reduces to 
tc ~ [S-y/vrGprp/cj] ^ . For the M33 nuclear parameters, this gives 

= 7 X 10^ f -V' ( yr. (9) 

V2 X IO^Mq pc-3/ V21kms-V ^' 

In other words, even for an intermediate-aged nucleus, the probability that a solar-type star will 
have experienced a collision at the center of M33 approaches unity. 

If we assume that the present turnoff of the parent old population in M33 is at ~ IMq, then 
the typical IMq blue-straggler collision product will resemble an A5 V star. Explaining the larger 
B — R versus V — I color gradient requires only that stars added to the center of the nucleus be of 
spectral type earlier than F5 or so. For an A5 star with B — Rk, 0.3, and V — I 0.2, a fractional 
contribution to the total nuclear light of ~ 30% is sufficient to explain both central color gradients; 
a ~ 20% contribution of AO stars or a ~ 40% contribution of FO stars would also work. 

Due to the low mass-to-light ratios of A- type stars, the implied mass contributions are 
modest and are well within the predicted mass fraction of collision products. For an assumed 
M/Ly ~ 2.5 for the old nucleus population (see below) and an A5 luminosity of I6L0 at V, the 
implied mass- fraction of blue stragglers is only ~ 1 x 10~^. Allowing for the finite lifetime of an 
A5 star, t ~ 1.5 x 10^ yr, compared to the nominal 7 x 10^ yr collision time still implies that 
the blue population required is only ~ 6% of that predicted by the naive application of equation 
(^. Explaining the blue population with somewhat more massive and hotter collision products, 
such as AO stars, would imply a still more modest collision efficiency. Even a population of blue 
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stragglers dominated by FO stars still fits within the naive predictions of the collision rate by an 
order of magnitude. 

The sharp radial drop in stellar density in the outer nucleus means that collisions will be 
effective only at the very center. The collision rate calculated above is for r < 0''03; at r O'.'l, 
the rate is lower by an order of magnitude. The strong dependence of collision rate with radius 
might explain both the scale and onset of the central nuclear color gradients. The typical stellar 
orbit in a potential similar to that of the nucleus is likely to have a factor of three ratio between 
its pericenter and apocenter distances (assuming an isotropic velocity ellipsoid). Any collision 
products generated within r < 0^'03 are likely to have significant diffusion out to r ~ O^'l but have 
little presence at still larger radii, in qualitative agreement with the color gradients observed. 

A qualitative argument suggests that the above collision rate may indeed be somewhat 
overestimated. If the M33 nucleus has undergone relaxation and collapse as suggested by 
Hernquist, Hut, &: Kormendy (1991), then significant mass segregation may have taken place at its 
center, suggesting that the modest velocity dispersion measured by Kormendy & McClure (1993) 
may be more representative of the star light rather than mass. The massive blue stragglers should 
be in equi-partition with the older stars, suggesting that the true velocity dispersion, and hence 
collision time scale, is higher. 

A bias of this sort may in fact be suggested by the low M/Ly ~ 0.4 calculated by Kormendy 
& McClure (1993). Models presented by Worthey (1994) show that the old nucleus population will 
most likely have M/Ly > 2; a conclusion that appears robust over a wide range of metallicities 
provided that the bulk of the population is at least 5 Gyr old, an assumption supported by 
the spectral synthesis work of Schmidt, Bica, & Alloin (1990) discussed below. The low M/L 
blue-straggler population will dilute this of course, but not by the factor of five required; the 
above model with 30% of the light in A5 stragglers, for example, has M/Ly = 1.8. To obtain an 
apparent M/Ly = 0.4 would require measuring a a that is ~ 2x too low, not unreasonable if 
much of the light comes from A stars with masses ~ 4x that of the typical star. Regardless, such 
arguments imply that tc at best is only ~ 2x longer than a naive application of equation ^ might 
suggest. This does not change the conclusion that collisions are important, as the calculation still 
has a factor of ~ 5 to spare. 

The possibility that collisions are responsible for the young population in M33 contrasts with 
the traditional conclusion that the M33 nucleus is a site of repeated episodes of star formation 
since its initial formation (see van den Bergh 1991). Given the inherent freedom in appealing to 
episodic star formation, we cannot argue that the collision scenario is by any means favored, and 
indeed it may have difficulty accounting for the UV spectral energy distribution of the nucleus. 
Schmidt, Bica, & Alloin (1990) use spectral synthesis to conclude that a ~ 10^ yr population 
contributes ~ 30% of the nuclear light, consistent with the collision model, but they also argue 
that a significant contribution from a very young 10'' yr component is also present. Massey et al. 
(1996) also infer a hot young component in the nucleus from far-UV photometry. Under the pure 
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collision scenario, such hot components would have to be stars produced in second-generation 

mergers of the initial blue stragglers, or perhaps stars formed from gas left over from mass-loss in 
the collisions. Such effects can only be evaluated by full models of the nucleus rather than the 
simple arguments presented here. We only argue here merely point out that, unless the collision 
time has been severely underestimated, collisions should have readily observable effects on the 
gross properties of the M33 nucleus. 

4. Toward a Unified Picture 

We conclude with an attempt to place the Local Group nuclei into a common evolutionary 
scenario. As the nearest examples of nuclei in giant galaxies, they may provide a view of 
phenomena in more distant nuclei, including AGNs. For example, the steep density profiles of M31 
and M32 place them in the "power-law" category of central density profiles according to Lauer 
et al. (1995). Understanding the origin of M31 and M32 may thus illuminate the origin of all 
power-law profiles, which typify small ellipticals and bulges generally; "power-law" galaxies may 
also serve as the building blocks of the shallower core-type profiles of massive ellipticals (Faber et 
al. 1997). 

The Introduction strongly emphasized the distinct properties of these three nuclei; in 
actual fact their inner density profiles prove to be broadly similar, with central densities 
~ W^'^^^'^ Mq pc~^ at 0.1 pc and radial slopes p ~ |.-i-5±o.5_ However, central dispersions differ 
enormously, with M31 and M32 in the range 200-800 km/s while M33 has a much lower cr = 21 
km s~^. This difference of course reflects the putative presence of black holes; M31 and M32 have 
estimated black hole masses in the range 3 X 1O^~^M0, while a black hole is undetectable in M33 
at the level < 2 x lO^M©. 

This tentative division into two types according to black hole content is heightened by adding 
the fourth large galaxy in the Local Group, the Milky Way. The Milky Way nucleus combines the 
properties of M32ibl in interesting ways: it has a stellar density profile that is rather similar to the 
other three, with a peak central density at 0.1 pc ~ 2 x IO^Mq pc~^, falling off radially beyond 
0.3 pc as r~^'^ (Haller et al. 1996). Like M31 and M32, it has high central velocity dispersion 
(t(0.1 pc) ~ 200 km s~^ and an estimated central black hole mass similar to M32 of 2 x IO^Mq 
(Haller et al. 1996). However, it also has a young central cluster of blue stars (see review by 
Morris &; Serabyn 1996), and in this regard it more resembles M33. Significantly — unlike any 
of the others — the Milky Way is barred, which strongly perturbs the central gas dynamics and 
appears to create the conditions for a steady inflow of gas into the central ~ 5 pc at a rate of 
roughly 0.03-0.05 Mq yr~^ (Morris &: Serabyn 1996). Bar aside, the Milky Way seems to fit well 
with M31 and M32. Indeed, from the standpoint of central structure, the Milky Way looks like a 
barred spiral galaxy with a transplanted M32 spheroid. For purposes of comparison, we therefore 
lump it with M32 and M31. 
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As a group, it may be significant that the stellar collision time scales for all four Local 
Group galaxies hover at or below values that would have barely detectable dynamical effects over 
a Hubble time. This suggests that the systems have evolved via similar paths to this state, or 
will do so eventually. One such scenario has been outlined by Murphy, Cohn, &; Durisen (1991, 
MOD), who have modeled the evolution of compact star clusters in the vicinity of seed black 
holes, including the effects of black hole growth and stellar relaxation, stellar mass loss, and stellar 
collisions (but not mass ejection by supernovae). The inner density profiles and black hole masses 
of their Models 3 and 4 at an age of ~ 10 Gyr are qualitatively consistent with M31, M32, and 
the Milky Way. Interestingly, the initial central density profiles in all MCD models start out much 
lower and flatter within 1 pc than the evolved profiles — the dense, steep profiles seen at ~ 10 
Gyr are due to the combined effects of stellar relaxation, loss-cone depletion, and adiabatic pulling 
by the growing black hole. Early evolution of the stellar profile is rapid owing to the initially low 
stellar velocity dispersion, which induces short stellar relaxation and collision times. As the black 
hole grows, however, the stellar velocity dispersion rises and the stellar population gets "stiffer;" 
relaxation shuts down and collisions become less frequent. This is the same stabilizing by the 
black hole that was invoked by Lauer et al. (1992) to postpone rapid dynamical evolution in M32. 
The black hole continues to grow at late times via stellar mass loss and loss-cone capture, but only 
slowly, and the systems after 10 Gyr appear to be in near-steady state, like M31, M32, and the 
Milky Way. 

An important question is why the initial conditions of M33's nucleus set its evolution on a 
different path from the that of the others. A possibility is that the M33 galaxy lacks a massive 
spheroidal component while the other three galaxies have one. If initial central density is the 
key to subsequent evolution, as the MCD models suggest, then the presence of a surrounding 
spheroid, and consequent increase in stellar mass loss and/or merger- induced gaseous infall, may 
be needed to boost nuclear density to the required level to grow a black hole quickly. So far, all 
detected massive central black holes are in galaxies with spheroids (Richstone et al. 1998) save for 
one Seyfert nucleus in the bulgeless dwarf Sdm galaxy NGC 4395 (Filippenko, Ho & and Sargent 
1993). A major task for the Hubble Space Telescope is to search for black holes in bulgeless 
galaxies. 

The lack of a central black hole in the M33 nucleus leads to its profoundly short central 
relaxation time, t,. ~ 3 x 10^ yr, several orders of magnitude shorter than that for M32 or M31 
(see Table 1^. We have advocated the short central collision time in M33 as a possible explanation 
for its blue center; however, the lack of any apparent collision effects in M32, with a central 
collision time only ~ 3x longer serves as a strong caveat that this hypothesis as simply stated 
may be incomplete. The short relaxation time implies a strong likelihood that the M33 nucleus 
has core-collapsed. The collision probability increases strongly during collapse, suggesting that 
the blue-straggler population was formed during the collapse, rather than under any sort of 
quasi-static conditions. Indeed, the short relaxation time of the nucleus already implies that any 
naive calculations based on the presumption that its structure has not evolved over the age of its 
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stellar population are likely to be in error. In contrast, the long relaxation times in M31 and M32 
argue for little recent structural evolution, suggesting that the naive collision calculations may be 
closer to the truth in these two nuclei. 

The longer native collision time in M31 compared to M32 makes the presence of a central blue 
excess there but not in M32 a compelling problem. Because the blue region is resolved, it must be 
produced mainly by starlight. This is interesting because all conceivable types of hot stars — blue 
stragglers, hot horizontal branch stars, young upper main-sequence stars — are short-lived. Hence, 
understanding the origin of this blue flux sets constraints on current or very recent events in M31. 
Since the main-sequence stellar collision time is longer in M31 than in M32 (see j j3.l| and § |3.2D , we 
must appeal to some other difference between the two nuclei. We consider some possibilities in 
turn: 



1. The velocity dispersion near the black hole is ~ 3x larger at a given radius in M31 than 



in M32, so stellar collisions are more violent. As noted in §3.1, this works against the 
blue-straggler hypothesis for the blue excess in M31, as blue stragglers are less likely to 
coalesce at these high speeds. The higher overall velocity dispersion in M31 thus does not 
favor the formation of hot stars, at least not by mergers. 

2. The estimated black hole mass in M31 is ~ lOx larger than that in M32, and the maximum 
Eddington radiation limit is larger by the same factor. If both black holes radiated 
occasionally at their respective limits, the black hole in M31 would evaporate the envelopes 
of giants to a larger radius, revealing more bare hot cores. According to Tout et al. (1989), a 
radiation bath of 10*^ K causes K giants to lose their envelopes. However, at the M31 black 
hole Eddington limit of lO^^L©, the effective temperature of the radiation field falls to this 
level at only 1000 AU, ~ 40x smaller than region of the blue excess (0.2 pc). This model 
would also require an ad hoc recent maximal burst of AGN activity in M31 within the last 
few 10^ yr. 



3. As noted in § |3.1| and §3.2, the central stellar population may be younger in M32 than in 
M31, (2 X 10^ yr vs. 5 x 10^ yr). By Equation (0), this age difference merely nullifies the 
difference in the collision times, which also differ by a similar factor; it does not explain why 
the fractional blue excess is so much larger in M31 than in M32. 

4. A cold disk of stars in elliptical orbits is required in the Tremaine (1995) model for the 
M31 double nucleus, while no such disk is seen in M32. Suppose that some fraction / of 
stars within 0.2 pc in M31 are orbiting in a cold disk with a^isk = 0.05vrot, which is typical 
of young disk stars in the Milky Way. At 0.1 pc, the random velocity dispersion in the 
disk, adisk, is ~ 60 km s~^, ~ 120 x colder than the isotropic velocity dispersion of 700 km 
s^^. The focusing enhancement of the collision cross section over the geometrical value is 
the Safronov number (^) given by = {V^^'^ /2adisk)'^ , where V^^"^ is the escape velocity 
from the stellar surface (600 km s~^ for the sun), and Udisk is the 1-d velocity dispersion of 
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the disk. In addition, the effective density of the disk stars is also increased by the factor 
~ (t/ a disk a-s they are compressed into the disk. Under these assumptions, the new colhsion 
time is reduced by the factor ~ 0~^, which for cJcfisfc = 60 km s^^, is a factor of 25. The 
new collision time is 2.4 x 10^ yr, and thus essentially all disk stars would suffer a low-speed 
collision, causing them to merge. Roughly a fraction / of the central light would be emitted 
with the effective temperature of blue stragglers, and thus / ~ 0.2 is needed to match 
the magnitude of the blue excess. Though promising, this model still does not reach the 
high color temperature of > 10^ K that seems to be needed to match the spectral-energy 
distribution of the blue excess. 

5. M31 possesses neutral gas like the Milky Way — unlike M32, which is barren. The region 
within 1 pc of the Milky Way black hole appears to be a site of episodic star formation, 
with a mass of stars formed in the last few million years of about IO^Mq (Morris & Serabyn 
1996). In M31, the color temperature of the blue excess indicates a main-sequence turnoff 
mass near AO and therefore an age of about 10^ yr if modeled as a single-burst population. 
The blue-excess V luminosity of 1.7 x lO^L© and resultant MjLy ~ 1 (Bruzual h Chariot 
1993) imply a mass of 10^ Mq, similar to the Milky Way central cluster. The radii of the 
clusters are also similar, a few tenths of a parsec. If allowed to age undisturbed, the young 
population centered on the Milky Way black hole would evidently look rather like the blue 
excess population in M31 in about 1 Gyr. This explanation begs the question of forming 
stars in such a region, where the random velocities exceed several 10^ km s"^ and the shear 
is intense, but the Milky Way appears to have solved that problem somehow, albeit with 
velocities about a factor of three smaller. It is also true that the center of M31 is relatively 
empty of neutral gas now, but an intermittent bar instability might bring gas to the center 
on periods of a few Gyr. It is intriguing that, in this model, the dynamical structure of M31 
might have looked quite different in the past. 

None of these proposals for explaining the blue excess population in M31 is without drawbacks. 
However, all illuminate the importance of carrying studies of galactic nuclei to the next stage, 
where all significant processes are modeled and related nuclei are compared and contrasted. We 
hope that the emerging body of high-resolution data on Local Group galaxies will stimulate 
modelers to match in detail the evolutionary history of these nearby Rosetta Stone objects. 

We wish to thank Doug Lin as well as our many collaborators on the Nuker team for 
useful discussions. This research was conducted by the WFPC-1 Investigation Definition Team, 
supported in part by NASA Grant No. NAS5-1661. 
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6 X 10 


1.6 X 


lO'' 









Table 2. 


M31 Nuclear Colors 




Color 


PI 


P2 Anti-Pl 


1" < r < 5" 


V-I 


1.41 ±0.01 


1.28 ±0.03 1.38 ±0.01 


1.34 ±0.01 


F300 - V 


2.29 ± 0.02 


1.29 ±0.03 2.18 ±0.02 


2.39 ± 0.02 



Note. — The PI and anti-Pl colors were measured in 0''4 x 0''4 
boxes centered on PI and displaced 0''4 from P2 in the anti-Pl 

direction, respectively. The P2 colors were measured in a smaller 
0''14 x 0('14 box to isolate the UV-bright portion of the nucleus. 
The last column represents the bulge outside the nucleus; the radial 
range is distance from P2. 
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Table 3. M31 Deconvolved Surface Photometry 





a 










i^300 




V 




V — 


1 






PA 




e 









0114 


13 


235 zb 


050 





763 zb 





070 


1 


169 zb 





070 

















0227 


13 


340 ± 


004 





658 zb 





004 


1 


240 zb 





005 


60 


9 zb 8.4 





254 ± 





037 





0455 


13 


529 ± 


007 





655 zb 





009 


1 


357 zb 





009 


60 


9 zb 3.1 





320 ± 





017 





0683 


13 


603 ± 


020 





921 zb 





023 


1 


383 zb 





025 


60 


9 zb 10.4 





484 ± 





087 





0910 


13 


621 ± 


018 


1 


143 zb 





025 


1 


377 zb 





025 


64 


6 zb 25.1 





613 ± 





260 





1137 


13 


643 ± 


014 


1 


246 zb 





026 


1 


371 zb 





020 


60 


4 zb 10.8 





671 ± 





126 





1365 


13 


683 ± 


012 


1 


308 zb 





040 


1 


374 zb 





016 


57 


8 zb 4.2 





671 ± 





049 





1592 


13 


709 ± 


009 


1 


356 zb 





025 


1 


374 zb 





013 


56 


8 ±4.1 





671 ± 





048 





1820 


13 


733 ± 


008 


1 


397 zb 





030 


1 


394 zb 





013 


54 


1 zb 3.5 





653 ± 





040 





2048 


13 


748 ± 


007 


1 


484 zb 





023 


1 


372 zb 





010 


48 


7zb4.5 





634 ± 





049 





2275 


13 


790 ± 


006 


1 


533 zb 





019 


1 


404 zb 





008 


46 


8 zb 1.3 





588 ± 





013 





2828 


13 


794 ± 


007 


1 


716 zb 





022 


1 


363 zb 





009 


46 


7 zb 3.4 





578 ± 





022 





3327 


13 


812 zb 


007 


1 


833 zb 





017 


1 


341 zb 





009 


56 


7 zb 6.6 





577 ± 





014 





3915 


13 


924 ± 


007 


1 


857 zb 





020 


1 


392 zb 





009 


53 


9 zb 4.9 





500 ± 





017 





4605 


13 


962 ± 


006 


1 


870 zb 





018 


1 


384 zb 





008 


57 


3 zb 1.4 





520 ± 





027 





5418 


14 


014 it 


006 


1 


838 zb 





023 


1 


322 zb 





007 


62 


4 zb 0.5 





536 ± 





007 





6374 


14 


227 it 


006 


1 


766 zb 





021 


1 


354 zb 





007 


65 


zb 0.6 





402 ± 





008 





7499 


14 


330 ± 


005 


1 


826 zb 





018 


1 


338 zb 





006 


62 


7 zb 0.5 





410 ± 





006 





8822 


14 


497 zb 


005 


1 


914 zb 





016 


1 


380 zb 





006 


62 


7 zb 0.4 





405 ± 





005 


1 


0379 


14 


695 ± 


005 


1 


934 zb 





014 


1 


296 zb 





006 


62 


4 zb 0.4 





361 ± 





005 


1 


2211 


14 


893 zb 


004 


1 


886 zb 





Oil 


1 


356 zb 





005 


60 


5 zb 0.4 





332 ± 





004 


1 


4366 


15 


108 zb 


004 


1 


847 zb 





013 


1 


337 zb 





005 


55 


1 zb 0.5 





270 ± 





004 


1 


6901 


15 


212 zb 


003 


1 


897 zb 





012 


1 


324 zb 





004 


58 


zb 0.4 





308 ± 





004 


1 


9883 


15 


356 zb 


002 


1 


976 zb 





009 


1 


334 zb 





003 


54 


7 zb 0.4 





294 ± 





003 


2 


0759 


15 


429 zbO 


002 


1 


970 zb 





Oil 


1 


386 zb 





003 


55 


3zb0.5 





223 ± 





003 


2 


4422 


15 


523 ±0 


001 


1 


954 ± 





009 


1 


372 ± 





002 


58 


2 ±0.3 





252 ± 





003 


2 


8732 


15 


653 ±0 


001 


2 


042 zb 





009 


1 


336 ± 





001 


55 


6 ±0.3 





219 ± 





003 


3 


3802 


15 


768 zbO 


001 


1 


981 zb 





009 


1 


344 zb 





001 


55 


0±0.6 





137 ± 





003 


3 


9767 


15 


861 zbO 


001 


2 


026 zb 





007 


1 


348 zb 





001 


48 


7 ±0.6 





135 ± 





002 


4 


6785 


15 


951 zbO 


001 


2 


024 zb 





006 


1 


338 zb 





001 


58 


8 ±0.5 





133 ± 





002 


5 


5041 


16 


026 zbO 


001 


















49 


6 ±0.7 





107 ± 





002 


6 


4754 


16 


092 ±0 


001 


















47 


3 ±0.4 





154 ± 





002 


7 


6181 


16 


151 ±0 


001 


















47 


8 ±0.3 





182 ± 





002 


8 


9624 


16 


219 zbO 


001 


















42 


2 ±0.2 





203 ± 





001 


10 


5440 


16 


306 zbO 


000 


















40 


1 ±0.2 





170 ± 





001 



Note. — V-band deconvolved surface brightness profile is given as a function of 
isophote semimajor axis length in arcseconds. The area dominated by PI has been 
excluded from the fit. 
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Table 4. M32 Deconvolved Surface Photometry 





a 












T 7" T 

V — I 




PA 








e 




100 X ^4 


0, 


.0161 


10, 


.959 ± 


0, 


.070 


1, 


174 ± 090 




















0, 


.0455 


11, 


.228 ± 


0. 


.008 


1, 


216 ± 009 


160. 


.3 ± 


10.5 


0. 


.146 ± 0, 


.027 


-0, 


.7± 1.5 


0, 


.0910 


11, 


.642 ± 


0, 


.004 


1, 


247 ± 006 


160. 


.3 ± 


1. 


.1 


0. 


.230 ± 0, 


.004 


-2, 


.9 ±0.8 


0, 


.1365 


11, 


.891 ± 


0, 


.007 


1, 


239 ± 017 


160. 


.3 ± 


1. 


.8 


0. 


,285 ± 0, 


.009 


0, 


.1± 1.2 


0, 


.1820 


12, 


.094 ± 


0, 


,008 


1, 


.255 ± 0.012 


160, 


,3 ± 


1, 


,7 


0. 


,286 ± 0, 


.009 


0, 


.0± 1.1 


0, 


.2275 


12, 


.235 ± 


0, 


,008 


1, 


.220 ± 0.013 


160, 


.3 ± 


2, 


.1 


0. 


,285 ± 0, 


.010 


2, 


.4± 1.1 


0, 


.2730 


12, 


.425 ± 


0, 


.007 


1, 


.223 ± 0.013 


160, 


.3 ± 


1, 


.3 


0. 


,259 ± 0, 


.006 


-0, 


.2 ±0.9 


0. 


.3185 


12, 


.568 ± 


0. 


.008 


1, 


.230 ± 0.013 


159. 


.8 ± 


1. 


.8 


0. 


,247 ± 0, 


.008 


-0, 


.1± 1.0 


0. 


.3640 


12, 


.686 ± 


0. 


.010 


1 


241 ± 018 


159. 


.7± 


2. 


.3 


0. 


,243 ± 0, 


.010 


-0, 


.7± 1.1 


0. 


.4095 


12, 


.787 ± 


0. 


.007 


1 


255 ± 014 


159. 


.7± 


1. 


.8 


0. 


,247 ± 0, 


.008 


-1, 


.5 ±0.8 


0, 


.4550 


12, 


.898 ± 


0, 


,008 


1, 


.246 ± 0.014 


159, 


.7± 


1, 


.7 


0. 


,247 ± 0, 


.007 


1, 


.9 ±0.9 


0, 


.5005 


12, 


.979 ± 


0, 


,008 


1, 


.225 ± 0.014 


159, 


.7 ± 


2, 


.0 


0. 


,255 ± 0, 


.009 


2, 


.9 ±0.8 


0. 


.5460 


13, 


.076 ± 


0, 


,008 


1, 


.251 ± 0.014 


159, 


.5 ± 


1, 


.5 


0. 


,259 ± 0, 


.007 


-0, 


.6 ±0.8 


0. 


.5915 


13, 


.164 ± 


0. 


.008 


1, 


.263 ± 0.012 


159. 


.7± 


1. 


.4 


0. 


,259 ± 0, 


.006 


-0, 


.5 ±0.7 


0. 


.6370 


13, 


.250 ± 


0. 


.007 


1, 


.264 ± 0.011 


160. 


.0 ± 


1. 


.3 


0. 


,259 ± 0, 


.006 


-1, 


.4 ±0.7 


0. 


.6825 


13, 


.354 ± 


0, 


.004 


1 


258 ± 008 


160. 


.0 ± 


0. 


.5 


0. 


,259 ± 0, 


.002 


-0, 


.6 ± 0.4 


0, 


.7280 


13, 


.426 ± 


0, 


.006 


1, 


250 ± 010 


160, 


.0 ± 


1, 


,1 


0. 


,262 ± 0, 


.005 


-1, 


.4 ±0.5 


0, 


.7829 


13, 


.509 ± 


0, 


,011 


1, 


.217 ± 0.014 


160, 


.3 ± 


1, 


.3 


0. 


,264 ± 0, 


.009 


0, 


.6 ±0.4 


0. 


.9211 


13, 


.746 ± 


0, 


,010 


1, 


.223 ± 0.013 


159, 


.1 ± 


1, 


.1 


0. 


,281 ± 0, 


.008 


-1, 


.2 ±0.7 


1, 


.0836 


13, 


.951 ± 


0. 


.010 


1, 


.248 ± 0.013 


157. 


.7± 


0. 


.9 


0. 


,290 ± 0, 


.008 


-0, 


.2 ±0.3 


1. 


.2748 


14, 


.203 ± 


0. 


.010 


1 


225 ± 013 


159. 


.1 ± 


0. 


.9 


0. 


,288 ± 0, 


.007 


-0, 


.1±0.5 


1, 


.4998 


14, 


.450 ± 


0, 


.009 


1, 


243 ± 012 


158, 


.9 ± 


0, 


.8 


0. 


,282 ± 0, 


.007 


-0, 


.3 ±0.1 


1, 


,7645 


14, 


.684 ± 


0, 


,009 


1, 


.239 ± 0.011 


160, 


.5 ± 


0, 


.8 


0. 


,293 ± 0, 


.007 


0, 


.1±0.2 


2, 


.0759 


14, 


.906 ± 


0, 


,008 


1, 


.234 ±0.011 


159, 


,0± 


0, 


,7 


0. 


,301 ± 0, 


.007 


-0, 


.1±0.2 


2. 


.4422 


15, 


.142 ± 


0, 


.008 


1, 


.223 ±0.010 


159, 


.6± 


0, 


.8 


0. 


,284 ± 0, 


.007 


-0, 


.1±0.3 


2. 


.8732 


15, 


.351 ± 


0. 


.007 


1, 


.225 ± 0.009 


157. 


.9± 


0. 


.7 


0. 


,293 ± 0, 


.006 


0, 


.0±0.2 


3. 


.3802 


15, 


.580 ± 


0. 


.007 


1, 


.233 ± 0.009 


159. 


.9± 


0. 


.7 


0. 


,276 ± 0, 


.006 


0, 


.1±0.3 


3, 


,9767 


15, 


.785 ± 


0, 


,006 


1, 


.244 ± 0.008 


159, 


,1 ± 


0, 


,7 


0. 


,279 ± 0, 


.006 


-0, 


.2 ±0.1 


4, 


.6785 


15, 


.997 ± 


0, 


,006 


1, 


.222 ± 0.008 


159, 


,2± 


0, 


.7 


0. 


,269 ± 0, 


.005 


0, 


.6 ±0.3 


5, 


,5041 


16, 


.194 ± 


0, 


,005 


1, 


.238 ± 0.007 


158, 


.2± 


0, 


.6 


0. 


,283 ±0, 


.005 


0, 


.2 ±0.2 


6, 


.4754 


16, 


.408 ± 


0, 


.005 


1, 


.219 ±0.006 


158, 


.6± 


0, 


.5 


0. 


,290 ± 0, 


.004 


0, 


.0±0.2 


7. 


.6181 


16, 


.672 ± 


0. 


.005 


1, 


.239 ± 0.006 


161. 


.6± 


0. 


.4 


0. 


,292 ± 0, 


.004 


0, 


.2 ±0.2 


8. 


.9624 


16, 


.898 ± 


0. 


.005 


1, 


.220 ± 0.006 


160. 


.2± 


0. 


.4 


0. 


,300 ± 0, 


.004 






10. 


.5440 


17, 


.175 ± 


0. 


.007 


1 


.243 ± 0.009 


158. 


.8± 


0. 


.5 


0. 


,284 ± 0, 


.005 







Note. — V-band deconvolved surface brightness profile is given as a function of 
isophote semimajor axis length in arcseconds. 
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Table 5. M33 Deconvolved Surface Photometry 





a 












V - 1 




PA 








e 




u. 


.UlDl 


lU. 


.oo / It 


u. 


n7n 


u. 


.yzo It u.uyu 
















n 

u. 




1 1 

11. 


/1 77 -1- 
.4 ( / 31 


u. 


not; 
.uzo 




u. 


879 -1- n n^i 

.0 / Z It U.Uol 


D, 


.2 ±38, 


n 
.u 


n 

u. 


nc;9 4- n 

UOZ It u. 


0^9 
.UoZ 


n 


.uy iu 


1 9 

IZ, 


.ODU It 


n 

u. 


090 

.uzy 


n 
u 


.000 It U.Uoo 


D, 


.2± 12, 


.z 


u. 


o'?a 4- 
.Uoo It u. 


ooa 

.UUo 


n 


. ioDO 


io, 


.uoy It 


n 

u. 


.UDZ 


n 
u 


as*? -1- n 071^ 

.000 It U.U / 


D, 


.2± 11, 


9 
.z 


n 
u. 


1 ^0 4- 
,ioy It u. 


.UOl 


n 


. ioZU 


1 

io, 


.DZo It 


n 
u. 


099 

.uzz 


n 
u 


0/10-1-0 O'^O 

.y4u It u.uoy 


1 71 
1(1, 


,0± 


4, 


.u 


u. 


11/14-0 
, 1 14 It U, 


ooa 

.UUo 


n 


.ZZ / 


1 /I 
14, 


.UoZ It 


n 

u. 


097 
.UZ ( 


1, 


OO/I -1- O"?/! 
.UU4 It U.Uo4 




y. 


.6± 


5, 


.0 


u. 


1 1 « 4- 
,110 It u. 


01 1 

.Ull 


u. 


.Z ( oU 


1 /I 

14, 


.oOo It 


u. 


098 


1, 


no^ -1- n fi/ifi 

.UUO It U.U4D 


1 7 
1 1 . 


.8± 


6, 


A 

.4 


u. 


1 1 « 4- n 
,110 It u. 


ni ^ 

.UlO 


u. 


.oloO 


1 A 

14. 


.040 It 


u. 


097 
.UZ ( 


1, 


.uoy It u.uo / 


91 

Zl. 


.7± 


5, 


/I 
.4 


u. 


, 1Z4 It U, 


ni 9 
.uiz 


n 


.OD'iU 


1 /I 
14, 


,y4z It 


n 

u. 


098 
,UZo 


1, 


DQ/I -1- 0/17 
.Uo4 It U.U4 ( 


91 

Zl, 


,7± 


4, 


A 

.4 


n 
u. 


1 '?7 4- 
, 10 < It u. 


01 

.UlU 


n 




10, 


,ZoD It 


n 

u. 


09f^ 
,UZO 


1, 


0*^9 -1- 0/1/1 
.UoZ It U.U44 


99 
ZZ, 


,0± 


3, 


/I 
.4 


n 
u. 


1/19 4-0 
, 14Z It U, 


ooa 

.UUo 


n 


.40DU 


10, 


,oy4 It 


n 

u. 


09Q 

,uzy 


1, 


OOfi 4- 0/17 
.UUD It U.U4 ( 


9Q 

zo. 


,0± 


6, 


9 
.Z 


u. 


1 ^0 4- 

,ioy It u. 


01 7 

.Ul ( 


n 
u. 


• OUUO 


10, 


.Duy It 


n 
u. 


.Uoo 


1, 


0/10 4- C\K'i 

.U4y It u.uoo 


99 

ZZ, 


.7± 


4, 


c 
.0 


u. 


1 4- 

,ioy It u. 


01 ■? 

.Ulo 


u. 




10. 


8QQ 4- 

.ooy It 


u. 


.uoo 


i. 


.UoU It u.uou 


9Q 
Zo. 


.0± 


4, 


n 
.u 


u. 


1 4- n 
,100 It u. 


ni 9 
.uiz 


0. 


.5915 


16. 


.020 ± 


0. 


.035 


0, 


.909 ± 0.051 


23. 


.0± 


4, 


.8 


0. 


,165 ±0, 


.014 


0, 


.6370 


16, 


,120 lb 


0, 


,033 


1, 


.003 lb 0.050 


22, 


,7± 


7, 


.0 


0. 


,186 ±0, 


.023 


0, 


.6825 


16, 


,199 lb 


0, 


,030 


1, 


.005 lb 0.045 


23, 


.0± 


7, 


.6 


0. 


,186 ±0, 


.025 


0, 


.7280 


16, 


.365 lb 


0, 


.030 


1, 


.084 lb 0.048 


22, 


.7± 


3, 


.4 


0. 


,180 ±0, 


.011 


0. 


.7735 


16, 


.473 lb 


0, 


.053 


1, 


.011 ±0.075 


20, 


.2± 


9, 


.0 


0. 


,180 ±0, 


.028 


0. 


.8190 


16. 


.634 ± 


0. 


.028 


0, 


.995 ± 0.056 


20. 


.2± 


3, 


.2 


0. 


,164 ±0, 


.009 


0. 


.8645 


16. 


.796 ± 


0. 


.031 


0, 


.993 ± 0.057 


20. 


.2± 


3, 


.3 


0. 


,164 ±0, 


.009 


0, 


,9100 


16, 


,969 lb 


0, 


,037 


0, 


.920 lb 0.070 


7, 


.4± 


3, 


.5 


0. 


,164 ±0, 


.010 


0, 


.9555 


17, 


,107 lb 


0, 


,036 


0, 


.939 lb 0.056 


7, 


.4± 


3, 


.8 


0. 


,177 ±0, 


.012 


1, 


.0010 


17, 


.130 lb 


0, 


.037 


0, 


.979 lb 0.059 


7, 


.4 ± 24, 


.3 


0. 


,177 ±0, 


.073 


1. 


.0465 


17. 


.037 ± 


0. 


.096 


0, 


.766 ±0.110 


179. 


.8 ± 16. 


.0 


0. 


,177 ±0, 


.049 


1. 


.0920 


17. 


.356 ± 


0. 


.031 


0, 


.930 ± 0.057 


0, 


.1± 


1, 


.4 


0. 


,147 ±0, 


.004 


1. 


.1375 


17. 


.430 ± 


0. 


.037 


0, 


.957 ± 0.064 


8, 


.7± 


6, 


.6 


0. 


,169 ±0, 


.020 



Note. — V-band deconvolved surface brightness profile is given as a 
function of isophote semimajor axis length in arcseconds. 
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Table 6. Time Scales at 0.1 pc 





P 


a 










tc 


Galaxy 


(Mq pC-3) 


(km s~^) 


e 




(yr) 




(yr) 


M31 "Blue P2" 


1 X 10^ 


720 


0.2 


7 


X 10^^ 


6 


X lO^o 


M31 "Outer P2" 


4 X 10^ 


720 


0.2 


2 


X 10^2 


2 


X 10^^ 


M32 


7 X 10^ 


240 


1.6 


3 


X 10^ 


2 


X 10^0 


M33 


2 X 10^ 


21 


220 


3 


X 10^ 


7 


X 10^ 



Note. — Relaxation times and stellar collision times are given for 
solar- type stars at 0.1 pc from the centers of M32ibl; these are calculated 
from the local mass densities, velocity dispersions, and Safronov numbers 
shown. The two M31 entries are for densities implied by the blue P2 
center directly, or extrapolated inwards from the outer P2 cusp (see 
text). 
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